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Abstract 

We consider the design of a non-local MonteCarlo algorithm for SU(3) lattice systems 
according to the idea of embedding the degrees of freedom corresponding to the center of the 
group Z(3). As a crucial ingredient to reach this goal, we present a practical implementation 
of a cluster algorithm for Z{3) systems with general random pair interaction. 
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I. Introduction 



Developing non-local upgrade mechanisms for SU(N) lattice models proved to be a rather 
hard task, because of the impossibility to carry on a straightforward generalization of Wolff 's 
embedding idea.[|[ There are apparently no such algorithms available for SU(N) lattice 
systems, at present. The relatively small effort devoted to this goal has probably to do with 
the fact that the dynamics of gauge fields takes only a small part of the total computer time, 
and the feeling is that there is little to be gained in optimizing it. Still, large statistical 
errors on the gluon dynamics may result when large lattices will be used to explore more 
deeply the scaling region. We take the point of view that sometime it will pay to have an 
efficient algorithm for the pure gauge dynamics and we want to argue in this paper that 
a strategy similar to Brower-Tamayo's || embedding may possibly work, as suggested by 
Wolff ||]. We present some preliminary steps to be taken to develop such an algorithm, even 
if their true value will be appreciated (hopefully) only after these ideas will be implemented 
in a realistic calculation. In order to simplify the presentation, we explicitly cover the N=3 
two-dimensional sigma-model, but there seem to be no obstructions to generalize the idea to 
the case of four-dimensional QCD at finite temperature, which is the kind of real application 
we have in mind. According to ideas of Polyakov going back to the 70's, it is precisely the 
center of the group which governs the deconfining transition, hence the Z(3) degrees of 
freedom should be responsible for a critical slowing down near the critical temperature. 



In sect.Q. we describe the simple embedding we have in mind for SU(3); in sect. Ill 



derive the algorithm and we give a few details on its implementation while in sect. IV 
present some preliminary numerical results. 



we 
we 



II. Z(3) embedding 



Embedding Z{3) degrees of freedom into a SU(3) model, like the non-linear sigma model, 
can be accompished by assigning a cross section of the fibre bundle Z(3) — > SU(3) — ► 
SU{3)/Z{3). For instance such a section may be chosen as 

fCe z(3) 
V e SU(3) -^C-ul , w , s 

V ; S \U€ SU{3)/Z{3) 

|Arg{TrU}| < vr/3 (2.1) 

Since the bundle is non-trivial, the section is only defined up to a set of measure zero, but 
this is fine for our purposes. The action transforms as follows (here and in the following 
< ij > implies that i and j are nearest neighbours): 

S = E |Re{TrVtVj} 

<ij> 
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^RejCiCjlTrUjUj} 



<ij> 



and it appears as a Z{3) model in the variables £ with random couplings. In the gauge 
case, we would end up with a Z(3) lattice gauge theory with random coupling. Near the 
continuum limit the matrix WUj fluctuates around the identity, and the system is effectively 
ferromagnetic. 

To exploit this embedding to beat critical slowing down at the continuum limit one 
needs an efficient algorithm for the Z(3) degrees of freedom. While a cluster algorithm 
for Z(3) spin systems is well-known in the real-constant-coupling case (a special case of 
q-Potts model with q=3), the random coupling case requires some modification. A general 
scheme for local pair interaction has been devised by Niedermayer § and it covers the 
present situation. We are going to derive a cluster algorithm along a slightly different line 
of reasoning, but the outcome will be the same. Instead of concentrating on transition 
probabilities satisfying detailed balance, we shall introduce an equivalent ensemble with 
bond variables as new degrees of freedom, like one does in the Ising model. There are two 
ways to define the action in terms of bond variables and these correspond to two special 
cases in Niedermayer's approach. The resulting cluster dynamics is non trivial in the first 
case while clusters are statistically independent in the second case. 

III. The algorithm 

Let us set the notation: £ 6 Z(3) is a generic spin variable, its values ranging on the set 
(l,u!,0), where u = exp(27ri/3) . The action we consider is the following 



<ij> 

where i,j range over pairs of next-neighbours. 

The coupling constants are subject to a limitation which is not essential for the 
algorithm to work, but will simplify the presentation, namely we assume that 



for all a condition which is likely to be met near the continuum limit. Now let us 

write the partition function in the following way: 





Argyll < 2tt/3 



(3.2) 



Z = II Z H 

<ij> 
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Zij = exp(/3Re{^ij})(5^ iCj 

+ exp(/?Re {fijjw})^,^. 
+ exp(/3Re{fi ij w})(5 Cii(5Cj 
Ignoring a common factor FJ exp(/3Re {f2ij}), we can rewrite Z as follows: 

* II (<U. + H Wft^O + ^ H^C;) (3-3) 

where 

Wij(z) =exp(/3Re{( Z -l)n ij }) 

On each link ij let us define 

W<j =mm(W ij (u),W ij (u)) 

W> 3 = max (Wij-M, Wy(w)) (3.4) 



^ = ^ if ^ > = ^) (3.5) 
. u> otherwise 



and 



Under the condition above on the phase of O (Eq. |3.2|) , it turns out that W < < 1; as a 
consequence, if we rewrite Z^ as follows: 

Zii = <W, +W*(1- S Ci>Cj ) + (W> - W^)^ - (3.6) 
we can introduce bond variables as in Swendsen- Wang's algorithm, to get either 

Z = J2 II W 1 - "V/> + ^,0^ + (W> - W$) 5 CiA ,]\ (3.7) 

C,n <ij> L J 

or 

Z = Y, II {<W&^ (! - + <Wc, «»V - W$)] + 5 niJ , Q W<\ (3.8) 
C,n <ij> k J 

As a first option we then have the following updating algorithm: first step, the bonds 
are switched on with probability (1 — W^), provided Q = Q (this first step is identical 
to Swendsen- Wang's). The second step is given by the spin update, which is governed 
by the term S nijj i, forcing the same spin value on each connected cluster. These cluster 
spin values are chosen according to a distribution Z c i ust which takes into account cluster 
interaction which occurs at the boundaries. The effective action for the cluster dynamics 



can be extracted from the term involving S nijt o. Let C denote a generic cluster, Qc hs spin; 
then we have 

Z dust = II W((cXc>) 
C,C> 
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where 

iedC,jedC , . 

W((c, CcO= II ( w« + (w« - <Wc 

<ij> v 7 

The clusters' interaction dictated by this formula can be now treated as an ordinary lattice 
system, e.g. by a heat bath method. What we have realized here is a transformation 
from the original Z(3) lattice system to another Z(3) system indexed by clusters. The 
implementation which we have already tested (see next section) is suggested by Wolff's single 
cluster algorithm. The assumption underlying such a choice is that in some equilibrium 
regime the existing spin values on the boundary of the cluster can be used to evaluate the 
probability distribution even without growing the other clusters. 



As for the second option, starting from Eq.3.8 the bonds are switched on with probability 

'i-w< ifC* = g 
< if Q ~ "IjQ 

, otherwise 

There is no interaction between clusters (which may be read from the fact that the coefficient 
^n,o is C-i n dependent), but the interaction among spins within a cluster is non trivial. 
Instead of trying to thermalize the system on each cluster, we can just apply a global 
rotation chosen independently at random on each cluster. Niedermayer showed that this 
move satisfies detailed balance and ergodicity is obviously satisfied since clusters with just 
one site are possible. 



The setup outlined here would not work in case Eq.3.2 is violated somewhere. It is 
however straightforward to modify the formulae to account for this. 



IV. Preliminary numerical results 



We have implemented the first option (single cluster) to simulate a two-dimensional 
Z(3) model with random coupling. The coupling f2jj is uniformly distributed in the region 

(1 < Re{ft} < |) n (|Arg{fi}| < tt/3) . (4.1) 

In order to count sweeps in a fair way, we defined a sweep as consisting of an update 
on f2 followed by Nha s cluster updates, with Nhits adaptively chosen in such a way that 
Nhits x < N c > ~ L 2 , where < N c > is the average cluster size and L is the lattice size. In 
this way the computer time for a sweep is essentially independent on j3. The autocorrelation 
time Tj n fhas been measured for the total magnetization < £ > using the formula given in 
||. We report the correlation length £ for < Re {CiCj} > projected at vanishing transverse 
momentum, together with Ti^and the average linear size of the clusters. Data have been 
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taken on a 128 x 128 lattice on several runs of 5000 sweeps each (Tab.[l]) and on a 200 x 200 
lattice with 2500 sweeps only (Tab.^). 

The numerical data show a clear sign of critical behaviour around (3 ~ 0.615. At the same 
time the autocorrelation time redoes not increase more than linearly in £, at least up to 
(3 ~ .610; beyond this value £ becomes a substantial fraction of the lattice size and restarts 
to grow more rapidly. To determine the dynamic critical exponent for the algorithm one 
needs a very high statistics, and our data are too preliminary to get a reliable conclusion. 
The sharp rise in r^near the transition, if confirmed, would mean that the algorithm is not 
as efficient as we would have expected, but it may still be of some practical value - since 
the actual values are not dramatically high. One also remarks that going to bigger lattices 
tends to improve the performance. 



V. Conclusions 



We have examined the feasibility of a cluster algorithm designed to embed Z(3) into 
SU(3) in the spirit of Brower and Tamayo. Clearly most of the work is still to be done. 
We have to combine the present Z{3) algorithm with some local update mechanism for the 
SU(3)/Z(3) degrees of freedom. This should carefully deal with the coset condition 2A - 
a heat bath technique would probably be best suited. For the application to gauge theory 
one has to further modify the algorithm, but some progress was done recently on gauge 
Z(2) and U{1) systems, and also this problem will hopefully be overcome. 

What will the overhead be for such kind of algorithm, say in the study of the deconfining 
transition? Judging form recent progress in implementing cluster algorithms on massively 
parallel machines^] jjj, the price to be paid should not be too high. 
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Table 1: 128 x 128 



(3 




'Tint 




.560 


2.3 


1.0 


4 


.580 


3.0 


1.0 


5 


.600 


4.7 


1.2 


8 


.602 


5.2 


1.3 


8 


.604 


5.7 


1.3 


9 


.606 


6.3 


1.4 


10 


.608 


7.7 


1.9 


11 


.610 


12.0 


3.5 


16 


.612 


~25 


~22 


32 



Table 2: 200 x 200 



a 




Tint 




.600 


4.7 


1.0 


7 


.602 


4.9 


1.2 


8 


.604 


5.5 


1.3 


9 


.606 


6.2 


1.5 


10 


.608 


7.6 


1.6 


11 


.610 


9.3 


2.5 


14 


.612 


18.0 


12.0 


23 


.613 


~25 


«21 


33 


.614 


~32 


«40 


41 



